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A semiclassical Quantum Hydrodynamic model has been derived by taking the 
moments of the Wigner-Boltzmann equation. For the first time, the closure has 
been achieved by the use of the momentum shifted version of all order quantum 
corrected solution of the Wigner-Boltzmann equation and that has considerably ex¬ 
tended the applicability of the model towards the low temperature and high density 
limit. In this context, the importance of the correlation and exchange effects have 
been retained through the Kohn-Sham equation in the construction of the Wigner- 
Boltzmann equation. The validity of the approach is subject to the existence of the 
Taylor’s expansion of the associated Kohn-Sham potential. 
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I. INTRODUCTION 


The theory of many body system can be broadly classified in two categories. First one 
is the microscopic view where each member is treated individually. Their placement in 
phase space, their mutual interactions, their micro helds are brought under the focus. In 
the second method, the identity of an individual becomes unimportant and the particles 
are treated as a collection. The properties of the system are obtained by averaging over 
the micro properties of the particles and are reflected through some directly measurable 
quantities like pressure, temperature as a function of space and time. The microscopic 
nature is investigated by kinetic equations like the many particle Liouville’s equation. On 
the other hand, the macroscopic nature is revealed through the hydrodynamic equations 
obtained by taking the moments of the kinetic equations like the Vlasov equation. The 
choice of the method entirely lies on the nature of the problem under investigation. 

The formulation of quantum hydrodynamics attracted considerable attention in the 

past[l|. The hrst attempt in this direction is concerned with the quantum version of particle 

and momentum conservation equations which have been constructed from the Schrddinger 

equation by expressing the wave function in terms of a density dependent amplitude and 

a phase factor. In this method 2-41 the quantum effect is revealed through the Bohm po¬ 
ll 

tential which does not have any classical counterpart. Another approach}^ is related to 
the phase space analysis which is connected with the Wigner function and the moments of 
quantum kinetic equations. This route has been extensively explored to develop quantum 
hydrodynamic equations GHS] . This technique gives rise to an infinite hierarchy of moment 
equations where any member of the group is connected to the next order moment. There¬ 
fore, in order to maintain the consistency with the physical condition a closure technique 
has to be devised. In this context, Gardnerj9| performed the moment expansion of the 
Wigner-Boltzmann equation near thermal equilibrium and high temperature limit to derive 
the QHD model with a h? order quantum corrected stress tensor 9l-ll l| . In the next stage, 
QHD is obtained using the entropy extremization technique [l2, 13 1. Later this principle is 
applied to the system of identical particles with a newly dehned quantum entropy in terms of 
reduced density matrix where the indistinguishibilty factor is incorporated}^. Application 


of Grad’s moment method 


15[ |. in the construction of hydrodynamic model, has been also 


reported. In this article we have followed the approach of Gardner and are able to introduce 
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the quantum correction to all orders to extend the range of applicability of the concerned 
hydrodynamic model. In addition to that, the quantum exchange and correlation effects 
which are absolutely necessary in this context, have been taken care of through the Kohn- 
Sham potential which converts a many body interacting system to a non-interacting 
one-particle system with an effective potential. The conventional Wigner equation deals 
with the Schrddinger equation of wave mechanics. In order to incorporate the effect of 
correlation and exchange we have replaced the Schrddinger equation with the Kohn -Sham 
equation which is a nonlinear version of the former. To calculate the average values in the 
hrst three moment equations, we have used the all order semiclassical equilibrium solution 


of the Wigner-Moyal equation 17(]. Quantum Hydrodynamics has a wide 
bility in the context of semiconductor {q. 


astrophysical systems 


20 |, 


181] physics, thin metal dims 19(| and even in the 


range of applica- 


21| and in nuclear p 


lysics 


and employed for analyzing the thermistor 23(| theory and superfluidity 2^ . 


22l | . Similar models have been proposed 


II. DERIVATION OF THE HYDRODYNAMIC MODEL 


The Kohn-Sham equation is the Schrodinger equation of the hctitious non-interacting 
particles where the potential is replaced by the effective external potential (Kohn-Sham 
potential) to incorporate the correlation and exchange effects loj 


d'ljj I II 

^ + (t>kh1p 


where (pkh is the Kohn-Sham potential. 
The corresponding Wigner equation is 


dt 


m dx dx dp ^ 




( 1 ) 


( 2 ) 


where f{x,p,t) is the single particle quasi-distribution function. 


C, = l/{2f{2] + 1)! 

The above equation is written in a normalized form, where we have dehned the following 

t 

ly/mP 


normalized variables. 
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X 


X = 


P = 


I 


m 


where I and l\/mf3 are the length and time scale of the system and (3 is the Boltzman 
constant. The normalized eqnation is obtained as 


dt 


P 


_ d(j)kh df ^ 

f)rr fi'rt ^ 


d_l 

dx dx dp 


where 


i=i 


A = 


(3) 


lh^/3 

mP 


is the small expansion parameter. This eqnation is valid only if it is possible to develop the 


potential energy in a Taylor series 


25l |. To obtain the hydrodynamic eqnation, we need 


to take the moments of the Wigner eqnation. The first three moment eqnations are 

^ ^ }_d^ ^ Q 

dt m dx 


(4) 


d{p) ^ _ ^ dcpkh 

dt dx'm dx 


(5) 


+ _^/^\ = (p) 94>kh 

dt 2m dx 2m^ m dx 


( 6 ) 


To calcnlate the average qnantities, we have used 


containing all higher order quantum corrections 


izl 


he following momentum shifted solution 


f{x,p,t) = n ^2j 

j>0 


where, 


W2 


exp(aoi 


p 

2 


cioi — A -h -B 


A = 

B = - 


^ /■ A (f(t)kh 

8 J \ ^ Q dx^ 


d(f^kh 

dx 


dx 


-1 


d^(j), 


kh 


dx^ 


dx 


(7) 

( 8 ) 
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and 


flu — (1 + 


<P(t)kh .-}, 

6 dx"^ 


( 9 ) 


For j > 1, 


W2j = exp{U2j) 


where 

Uy = 'tA^’a^Ar 

i=0 V ^ 

Qij can be obtained from the following equations 


1 dOf^j i -|- 1 

2* dx 2* dx 


( 10 ) 


where 


bij 




in which, Di = coefficient of p2*+i Qf Hermite polynomial He 2 j+i{p) 
These equations are true for i=0 to i=j-l 


For i=j, 


( 11 ) 


1 da 




■ Q^2H1 


( 12 ) 


2* dx 

This is a hrst order equation and can be easily solved to obtain the value of ajj. Using this 
result, all the equations of this group can be successively solved to find out the remaining 


coefficients. This process has been illustrated in 


17|. 


Expanding the exponential factors of hF 2 , we can write 


flOl — 0A:h + V 

and 

-Any = -pV2 + 6p^ 

7] and 9 are containing all higher order terms obtained from ooi and an. 

Now, 

2 J 2?' 

f{x,p) = Aexp{-^ + (j)kh)exp{p + Op'^) ^ exp[A^^(^ a^j^)] 

^ j>l i=0 ^ 
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which can be written in the following form 


/2 


p 

f(x,p) ^ Aexp{—— + 4>kh)exp{'q)G 


(13) 


and 


QT. 


2n 


G = Y.-p 

n nl 


A2E,>iJELo--n,>inL 


j>ilh=oa»j i 




■p 


where exponential factors have been expanded as 


ex: 


p(0p'2) ^ ^ 


{ep‘ 


/2\n 


exp[A^^ 


2i 


n\ 


ij J 


(14) 


n and riij are positive integers, and A{x,t) is a slow varying function of x and t and 


p = p — mu 


u and p are the macroscopic fluid velocity and the molecular momentum respectively. 
Therefore, we can write 


< pp/m >= u < p > —P 


(15) 


where 


<p'p'/m>=—P (16) 

and 

< pp^/2m‘^ >= uW — uP + q (17) 

where P and W are stress tensor and energy density respectively. 

q =< p'p'"^/‘2w? > (18) 

is the heat flux term and carries important effects of higher moments. 


W =< p^/2m >= mnu‘^/2+ < p'‘^/2m > 


( 19 ) 
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Using eq. ffT^ we obtain 


where 


< p >= mnu 


n{x,t) = j fix,p)dp 

Performing the momentum integral with the help of the following identity, 

/ x‘^"'e~^dx = \/n- 

Jo 

we obtain 


( 20 ) 


( 21 ) 


„2n„-4 l)'[^2n+l 


2 ' 2 n+l 


n[x 


on rr. in-’ •?' 

, t) = AV^expiv) E ^ ^ ^(2n + 2 ^ ^ m., - 1)!(22) 


2E,>iELo--n,>inLon*i! ^ 


J>1i=0 


Similarly, 


an A 2 Z3>i ■?' ZLo n ■ 1 n-iZ , /m 

P = AV^exp{p)Y ,—^^— . -——- ^(2n + 2^^my + I)!! (23) 


2E,>iZLo--n,>inLo’^*.! ^ 


j>li=0 


We define 


c = E 


Qn ^2 Zj>i ] J2i^o njj ^ 


(2n + 2^^my - I)!! 


j>li=0 


^ ^ ^ + 2 E E . D! 




J>1i=0 


Therefore 


< p >= mnu 

P = nB/C 

W = mnu‘^/2 — nB I2C 
Inserting < p >, P and W in flT^ and (ITTD 

< pp/2m >= [mnu^ — nB/C) 


(24) 


(26) 


(26) 

(27) 

(28) 

(29) 







< pp^/2m^ >= u{mnu^/2 — 3nB/2C) + q (30) 

As pointed out by Gardner we have allowed a heat flux term q = —k'VT, even though 
our distribution function does not allow it. Eq.(jl])-(|6]) constitute the set of hydrodynamic 
equations where < p >,< pp/2m > and < pp^/2m? > are obtained from eq. (l26ll . (|2^ and 
fl30|l respectively. 


III. CONCLUSION 


In this article we have established the Quantum Hydrodynamic model that incorporates 
quantum corrections to all orders. To do this we have followed the procedure introduced by 
Gardner who has developed the 0(/i^) fluid model by taking the moments of the Wigner- 
Boltzman equation to obtain the particle, momentum and energy conservation equations. 
In this process, Gardner has used the local ’’momentum-shifted” version of the 0{h?) solu¬ 
tion of Wigner equation with slow space and time dependence of density, temperature and 
macroscopic velocity to calculate the associated average values with the expectation that 
the equations will be valid, as in classical hydrodynamics, to describe systems even not near 
thermal equilibrium condition. Here, we have introduced the ’’momentum-shifted” version 
solution containing all higher order quantum corrections. Gonsequently, all higher order 
quantum terms are introduced to the stress tensor and the energy density and the range of 
validity of the QHD model has been extended considerably towards the high density and 
low temperature limit. In this connection, it is important to include the exchange and cor¬ 
relation effects of the particles and that has been taken care of through the application of 
the Kohn-Sham equation which has replaced the Schrodinger equation in the construction 
of the Wigner equation. This should be emphasized here that as we consider the higher 
order quantum terms, the correction due to the exchange and correlation should become 
more and more important and that implies the necessity of considering those effects when 
we are interested in the Quantum Hydrodynamic model with all order correction terms. 
The condition of validity of this model is restricted to the existence of the proper Taylor 
expansion of the potential function of the Kohn-Sham equation. 

At the end, it should be mentioned that in this model the distribution used to calculate 
the average quantities is unable to produce the heat flux. One way to get around this 


problem is the application of the Ghapman-Enskog expansion 


26| of the kinetic equation 
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with the BGK type relaxation term at the right hand side. This problem will be taken up 
as our future work. 

Finally, in the context of Quantum hydrodynamic model, the importance of this work is 
the construction of a more complete quantum fluid model with an extended probing capacity 
for analyzing low temperature and high density systems of different branches of physics. 
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